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NUMERICAL CALCULATION OF POTENTIAL-ENERGY CURVES 


BY RYDBERG-KLEIN-REES METHOD 
"by Frank J. Zeleznik 
Levis Research Center 


SUMMARY 

A technique is presented for the numerical e-valuation of the integrals oc- 
curring In the Rydberg-Klein-Rees method of calculating potential energy curves* 


INTRODUCTION 


The potential -energy curves for the hound states of diatomic molecules can 
he obtained from spectroscopic constants by using the Rydberg-Klein-Rees method 
(refs. 1 to 3). Klein (ref. 2) expressed the turning points of motion r max 
and r in terms of two auxiliary functions f and g: 
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Both of the functions f and g depend parametrically on the potential 
energy U and an additional parameter K = J(J + l)h /2p where J is the 
rotational quantum number and \i is the reduced mass. The functions f 
and g are defined as 
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Here E(l,K) represents the vibrational -rotational energy levels , and 
I = h(v + l/2) with v representing the vibrational quantum number. From 
spectroscopic data, the energy levels are generally expressible in the form 
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where are the usual spectroscopic constants. The upper limits of inte- 

gration in equations (l) and (2) are obtained from E(l',K) - U = 0. 

Rees (ref. 3) evaluated the integrals of equations (l) and (2) analytical- 
ly for special cases -where E(l,K) was either quadratic or cubic in I. The 
results in the cubic case were not in a very convenient form for computation. 
The results in the quadratic ease, however, have been used as the basis for 
rather extensive calculations by Vanderslice and coworkers (refs. 4 and 5) and 
others (refs. 6 and 7). Molecules, whose data could not be adequately repre- 
sented over the entire range by a quadratic, were treated by piecewise fitting 
of quadratics to the energy levels. As pointed out by Weissman, Vanderslice, 
and Battino (ref. 8), this piecewise fitting can lead to errors. 

The difficulty in the numerical integration of equations (l) and (2) is 
caused by the fact that the denominator of the integrands has a zero at the 
upper limit of integration. Jarmain (ref. 9) -circumvented this problem in an 
approximate manner by fitting [U - E(l,K)] 1 / 2 to an expression of the form 
c(l* - l)~^. The constants c and d were evaluated by using two points very 
close to the upper limit of integration. Similarly, Se/Sk was fitted to a 
quadratic in I 1 - I. These approximations were then analytically integrated 
to evaluate the contributions to f and g from regions close to I T . More 
recently, Weissman, Vanderslice, and Battino (ref. 8) introduced a new integra- 
tion variable x = [U - E(l,K)V 2 in order to eliminate the singularity from 
the integrands. Although correct, this procedure produces some unnecessary nu- 
merical inconvenience since the integrands are available as expressions in I 
and not the integration variable x. Thus, in numerical evaluation of the in- 
tegrals one cannot use arbitrarily selected increments of x but rather must 
select increments of I and calculate increments of x as 

Ax = [U - E( I + AI^K)] 1 / 2 - [u - e( i,k) ] V 2 

The technique that will be described here does not require the fitting of 
the integrand near the upper limit, and further, it essentially retains the 
original integration variable. It is based on the fact that the singularities 
in the integrands of f and g can be easily removed by an integration by 
part s . 


ALTERNATE METHOD FOR NUMERICAL EVALUATION 

The quantity U - E that appears in the integrands of f and g is a 
polynomial of order p in the variable I with coefficients that depend 
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on K. By introducing the notations x = v + l/2 and K = j(j + l) , the 
following dimensionless form is obtained: 
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= P(x;k) = A|(k)x^ 
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where E Q is a constant with the dimensions of energy and where the coeffi- 
cients in the polynomial P have the explicit form 
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In this notation the upper limit of integration I 1 in equations (l) and (2) 
corresponds to hx-i (k) where x-»(k) is the smallest, real, positive root of 
P(XJK). 


If the existence of the improper integral f is assumed, then equa- 
tion (l) Implies that the polynomial P(x;k) has a zero of order one at x-^(k). 
If Kq is a particular value of K, then a change in scale z = x/x-j_(kq) can 
be made, and the polynomial P can be written as 
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where A ^(k) = x^(Kq)A^(k) and where R(zjk) is a polynomial in z of order 
p - 1. If K is chosen equal to Kq, then equation (4) takes on the simple 
form 
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where the subscript has been dropped from k q for convenience, and where now 

a z (k) = x i(' c ) A j( K ) 
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Comparing coefficients of like powers of z in equation (5) establishes the 
relation between and as 
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Also from equation (5), for 


z = 1 one can obtain 
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Combining this with equation (6) finally gives the result 

P 

B z (k) = - ^ Aj(ic) l = 0, 1, . * p - 1 ( 7 ) 

j=Z+l 

The numerator of the integrand of g is essentially Bp/S/c. This partial 
derivative is most conveniently evaluated by using the first equality in equa- 
tion (4). After the differentiation is performed, k is again chosen to be 
equal to K q . This gives, after dropping the subscript on Kq, 
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The integrals for f and g can be written in dimensionless form by 
using the Bohr radius aQ as a unit of length and Eq = h^/2qa^ as a unit of 
energy. Also, making a change in scale to the new variable z and considering 
the particular value K that is equal to Kq result in 
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where the subscript on /c Q has again been deleted for convenience. Substi- 
tuting equation (5) for P(z ;k) into equations (9) and (10) and performing an 
integration by parts give 



In equations (ll) and (12) the objective of removing the singularity of the 
integrands at the upper limit has been achieved. 

For the important case k = 0, the polynomial coefficients required to 
evaluate the integrands of equations (ll) and (12) have the relatively simple 
forms 
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No difficulties were encountered in the numerical application of equations 
(ll), (12), and (13) to the calculation of potential-energy curves for several 
diatomic molecules by using standard integration techniques. In fact, the re- 
sults of single and double precision calculations for the ground state of hy- 
drogen agreed to at least seven figures when using the spectroscopic constants 
of Weissman, Vanderslice, and Battino (ref. 8). 


Lewis Research Center 

National Aeronautics and Space Administration 
Cleveland, Ohio, September 18, 1964 


REFERENCES 

1. Rydberg, R. : Graphic Representation of Results of Band Spectroscopy. Zs. 

Phys. , vol. 73, nos. 5-6, Dec. 10, 1931, pp. 376-385. 

2. Klein, 0. : Calculation of Potential Curves for Diatomic Molecules with Help 

of Spectral Terms. Zs. Phys., vol. 76, nos. 3-4, May 31, 1932, pp. 226- 
235 . 


5 



3. Rees, A. L. G. : The Calculation of Potential-Energy Curves from Band- 

Spectroscopic Data. Proc. Phys. Soc. (London), vol. 59, Nov. 1947, pp. 
998-1008. 

4. Vanderslice, J. T. , Mason, E. A., Maisch, W. G. , and Lippincott, E. R. : 

Ground State of Hydrogen by the Rydberg-Kle in-Rees Method. Jour. Molecu- 

lar Spectroscopy, vol. 3, no. 1, Feb. 1959, pp. 17-29. (See also errata, 
vol. 5, no. 1, July 1960, p. 83.) 

5. Fallon, Robert J. , Vanderslice, Joseph T. , and Cloney, Robert D. : Potential 

Curves and Rotational Perturbations of CN. Jour. Chem. Phys., vol. 37, 
no. 5, Sept. 1, 1962, pp. 1097-1100. 

6. Verma, R. D. : Ultraviolet Resonance Spectrum of the Iodine Molecule. Jour. 

Chem. Phys., vol. 32, no. 3, Mar. 1960, pp. 738-749. 

7. Singh, N. L. , and Jain, D. C. : The Rydberg-Klein-Rees Method of Construct- 

ing the True Potential Energy Curves of Diatomic Molecules. Proc. Phys. 
Soc. (London), vol. 79, Feb. 1962, pp. 274-278. 

8. Weissman, Stanley, Vanderslice, Joseph T. , and Battino, Rubin: On the Re- 

calculation of the Potential Curves for the Ground States of I 2 and H 2 . 
Jour. Chem. Phys., vol. 39, no. 9, Nov. 1, 1963, pp. 2226-2228. 

9. Jarmain, W. R. : Transition Probabilities of Molecular Band Systems. XX. 

Tabulated Klein-Dunham Potential Energy Functions for Ten States of Co(4), 
0g(2), OH ( 2 ) , and Si0(2). Sci. Rep. 8, Univ. Western Ontario (Canada), 
July 1, 1961. 


6 


NASA -Langley, 1964 E-2606 


"The aeronautical and space activities of the United States shall be 
conducted so as to co?itribute ... to the expansion of human knowl- 
edge of phenomena in the atmosphere and space. The Administration 
shall provide for the widest practicable and appropriate dissemination 
of information concerning its activities and the results thereof 

— National Aeronautics and Space Act of 1958 


NASA SCIENTIFIC AND TECHNICAL PUBLICATIONS 


TECHNICAL REPORTS: Scientific and technical information considered 

important, complete, and a lasting contribution to existing knowledge. 

TECHNICAL NOTES: Information less broad in scope but nevertheless 

of importance as a contribution to existing knowledge. 

TECHNICAL MEMORANDUMS: Information receiving limited distri- 

bution because of preliminary data, security classification, or other reasons. 

CONTRACTOR REPORTS: Technical information generated in con- 

nection with a NASA contract or grant and released under NASA auspices. 

TECHNICAL TRANSLATIONS: Information published in a foreign 

language considered to merit NASA distribution in English. 

TECHNICAL REPRINTS: Information derived from NASA activities 

and initially published in the form of journal articles. 

SPECIAL PUBLICATIONS: Information derived from or of value to 

NASA activities but not necessarily reporting the results of individual 
NASA-programmed scientific efforts. Publications include conference 
proceedings, monographs, data compilations, handbooks, sourcebooks, 
and special bibliographies. 


Details on the availability of these publications may be obtained from: 

SCIENTIFIC AND TECHNICAL INFORMATION DIVISION 

NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 

Washington, D.C. 20546 



